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Abstract: In this paper, we obtain an explicit formula for the two-point correlation 
function for the solutions to the stochastic heat equation on M. The bounds for p-th 
moments proved in [3] are simplified. We validate the Feynman-Kac formula for the 
p-point correlation function of the solutions to this equation with measure-valued 
initial data. 
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1 Introduction 

Consider the following stochastic heat equation 

{ f d v \ 

“ 2 ^ j ^ xeR, t>0, 

“(0, •) =h(-) , 

where lU is a space-time white noise and p : M —)■ M is a globally Lipschitz function. The 
initial data p is a signed Borel measure, which we assume belongs to the set 

:= |signed Borel measures p, s.t. J e^'^'^^IpKdx) < -foo, for all a > o|. 
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In the above, we denote |/i| := //+ + /i_, where /r = /i+ — /i_ and are the two non-negative 
Borel measures with disjoint support that provide the Jordan decomposition of /i. The set 
can be equivalently characterized by the condition that 

{\fi\ * G^{t, ■)) (x) = / G^{t, X — y)\^\{dy) < +00 , for alH > 0 and x G M, (1.2) 
Jr 


where * denotes the convolution in the space variable and Gy{t,x) is the one-dimensional 
heat kernel function 

The initial condition m( 0, •) = p(-) is understood as hm 4 oM(^, •) = in the sense of 

distribution (we identify a measure as a distribution in the usual sense; see [4, Theorem 

1.7)). 

Denote 

Jo(f,x) ■= {y,*Gy{t,-)){x) = / Gy{t,x -y)y{dy). 

Jr 

Dehne the kernel function 


/C(f, x) = /C(f, x; z/, A) 


G|(f,x) 


\/47rz/f 


A'^ Ah I 

-|-e 41^ <|> 

2v 



(1.3) 


where <h(x) = f^^(27r) Some functions related to $(x) are the error functions 

erf(x) = Jq e~y^dy and erfc(x) = 1 — erf(x). Note that <I)(x) = erfc(—x/\/2)/2. 

When p{u) = Xu, the following moment formula is proved in [3] 

E[u{t,xy) = Jo(t,x)-h (Jq */C)(f,x), (1.4) 

where ‘V’ denotes the convolution in both space and time variables, that is, 

{Jl-klC){t,x) ■.= f ds f dy Jl{s,y)lC{t - s,x - y). (1.5) 

Jo Jr 

As for the two-point correlation function, dehne 

X{t, xi, r, xa; v, A) :=A2 [ dr [ dz [Jq (r, z) + (Jg (•, o) * /C(-, o; z/. A)) (r, z)] 

Jo Jr 

X Gy{t — r,Xi — z)Gy{T — r,X 2 — z). 


Then by [3, (2.26)], for all r > t > 0 and xi,X 2 G M, 

E [u{t, xi)u (r, xa)] = Mt, Xi)Jo (r, Xa) + X(t, Xi, r, Xa; z^. A) . (1.6) 
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The first goal of this note is to simplify the moment formulas (1.4) and (1.6) for the 
case t = T. Note that the terms (Jq * )C){t,x) and X(...) involve four and six integrals, 
respectively. We will reduce these integrals into only two integrals: two convolutions of the 
initial data with respect to a kernel function. 


It is well known that if the initial data is a function, then the moments of the solution to 
(1.1) with p{u) = \u admit a Feynman-Kac representation; see [7]. Suppose that fi{dx) = 
Uo{x)dx where Uq is a bounded measurable function. This representation says that for all 
Xj G M, i = 1,..., n. 


E 


Y[u{t,Xi) 


2=1 


= E^ 


n 

2=1 


Mo [Xi 


Bl) 


exp 


E {Bi - Bi) ds 


l<i<j<n 


, (1-7) 


where {Bl,t > 0}, i = 1,... ,n, are i.i.d. standard Brownian motions on M, and the expec¬ 
tation is with respect to all these Brownian motions. 

On one hand, direct evaluation of this expectation when n = 2 is not easy since it involves 
the joint law of a standard Brownian motion Bt and its local time at an arbitrary level 
a G M. To the best of our knowledge, we are not aware of any references for this joint 
law expect the case where a = 0. We will derive this joint distribution and then give an 
alternative and more probabilistic proof of the formula for the two-point correlation function. 

On the other hand, when the initial condition is a measure the meaning of (1.7) is not 
clear. Another aim of this paper is to make sense of (1.7) for initial data in A4 h(E) using 
Malliavin calculus. More precisely, we show that nr=i + Bl) belongs to the Meyer- 
Watanabe space of Wiener distributions for any p > 1 and a > n{l — 1/p), and the 

exponential factor in (1.7) 

y„:=exp(A^ f Sr,-,. {Bi - Bi) dPj (1.8) 


belongs to D“’^ for any p > 1 and a < 1. Then, we can choosing p such that n(l — 1/p) < 
a < 1/2, we can write 


E 


Y\_u{t,Xi) 


i=l 


JJuo {xi + Bl) .Y,. 


(1.9) 


I 2 = 1 


where (•, •) denotes the duality between © and if 1/p -|- 1/g = 1. 

When p{u) in (1.1) is nonlinear but satishes the global Lipschitz condition, the explicit 
formula for the moment of the solution is impossible. We obtain an upper bound for the 
p-moment of the solution and a lower bound for the second moment. The idea is to compare 
them with the ones in the linear case. 


We hrst state our main results in Section 2. The result for the two-point correlation 
function. Theorem 2.1, is proved in Section 3. In Section 4, the joint law of {Bf, L“), Theorem 
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2.8, is proved. In Section 5, we give the alternative proof of our two-point correlation formula 
for function-valued initial data. Finally, by proving Theorems 2.10 and 2.11 in Sections 6.1 
and 6.2, respectively, we make sense of (1.7) for measure-valued initial data. 

Throughout of the paper, H-Hp denotes L^iVL) norm. 


2 Main Results 

2.1 Formulas for two-point correlation function 

Theorem 2.1. Suppose that pi G A4h(M)- If p{u) = \u, then for all t > 0 and xi,X 2 G M, 
E[u(t,Xi)u(t,X 2 )] = // pi{dZi)fi{dZ 2 ) JC*(t,Xi - Zi,X 2 - Z 2 ,Xi - X 2 ; \), (2.1) 


or 


E [u{t, xi)u{t, X 2 )] = Joit, xi)Jo{t, X 2 ) 

+ 11 pi{dZi)p{dZ2)K){t,Xi- Zi,X2- Z2,Xi-X2-,\), (2.2) 


where 


K){t,zi,Z 2 ,y,\) 2 [^^^- 2 (| 2 /| + \y-{zi - 22 ) 1 )]^ 


X <F 


i\y\ + \y- (^1 -^2)1) 

\/2vt 


(2.3) 


and 


/C*(t, zi, Z 2 , y; A) =G^{t, zi)G^{t, Z 2 ) + IC^it, zi, Z 2 , y; A) 

Zl+Z2^ G2y{t,Zi-Z2) 


— Gjy/2 

+ ^exp [XH-2{\y\ + ||/- (zi - ^ 2 )!)]^ 

- (bl + \y- (^1 - ^2)1)' 


(2.4) 


X $ 


\/2vt 


In particular, 


u{t,x)\\\ = // pi{dzi)pi{dz 2 ) lC{t,x - Zi,X - Z 2 ,d-.,\). 


(2.5) 
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In the following, we will use the convention that for any pair {wi,W2) of variables, 

w ■.= {wi + W2 )and Aw := W2 — Wi. (2.6) 

Remark 2.2. Formulas (2.1) and (2.2) are in the convolution form. One can also write 
them in the following inner product form: 


E[u{t,xi)u{t,X2)] = iJ,{dzi)fi{dz2) K* {t,xi,X2,zi,Z2 ;X), (2.7) 


or 


E[u{t,Xi)u{t,X2)] = Jo{t,Xi)Jo{t,X2) + // IJ.{dZi)lJ,{dZ2)K^t,Xi,X2,Zi,Z2-,X), (2.8) 


where 


u w f A^(|Aa;| + |A2;|) XH 

K' {t, xi,X2, zi, Z2; X) =—G^/2{t, X- z) exp I - — -h — 

X«. ( 1 )^ - . 

V \/2id ) 


(2,9) 


and 


K* (f, xi, X2, ^1, Z2) =Gi,(t, xi - zi)Gi,(t, X2 - Z2) + {t, xi, X2, zi, Z2; X) 


=G^/2 (t,x- z) 


G2u{t, Ax - Az) + — exp 


A2(|Ax| + |A^ ^ XH 


2z/ 


|Aa:| + |Az| \ 
V \/2z/ ^/2lyt ) 


Example 2.3 (Delta initial data). When fi = 60, then 


and 


\u{t^x )\\2 = IC*{t, X, X, 0; A) = A ^ /C(f, x). 


E [u{t, xi)M(f, X 2 )] =IC*{t, xi, X 2 , xi - X 2 ; A) 

=GHt, xi)GHt, X 2 ) + )C\t, xi, X2, xi - X 2 ; A) 

=G^(t,xi)G^(t,X2) + X— (t, —— 


X exp 


2l> ^ y 

X‘^{XH — 2|xi — X 2 I) 


$ 


X^t — |xi — X2I 
yj 2 vt 


4 . 


( 2 . 10 ) 


This recovers the results in [3, Corollary 2.8]. 
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Remark 2.4. Note that the function m 2 {t,xi,X 2 ) = K*{t,Xi,X2,0,0), or equivalently 
m 2 {t,xi,X 2 ) = ]C*{t,xi,X 2 ,Xi — X 2 ), solves the following parabolic equation (see e.g., [2, 
Theorem 3.2 on p. 46]) 

{ d 

—m2(t, xi, X 2 ) = H2{v, A) m2(f, xi, X 2 ), t> 0, xi, a;2 G M, 

7712(0, xi, 0:2) = ( 5 o(xi)(fo(a; 2 ), 
where the operator 

is the 2-particle Schrodinger operator (see [2]). 

Example 2.5 (Lebesgue’s initial measure). When /i(dx) = dx, then from (2.2) or (2.8), 


E [n(f, xi)M(f, X2)] = 1 + // dzidz2)C\t,zi,Z2,xi - X2;X) 

y 2 /• / y 2 \ 

= 1 + — y d^exp f — [A^t - 2 (|xi - X2I + |xi - Xg - z\)]\ 
'XH - (|xi - X2I + |xi - X2 - 2)1)' 


X $ 


\/2vt 


= 2e 41. 


- 2 A 2 ni-X 2 | / — |xi — X 2 I 


■<h 


\/2vt 


+ 2 $ 


|Xi - X 2 I 

yj2vt 


and in particular. 


|M(f,x )||2 = 2e 41. 




1, (2.11) 


( 2 , 12 ) 


These two formulas (2.11) and (2.12) recover the results in [3, Corollary 2.5]. Note that the 
equality in (2.11) can be established through integration by parts. 

Corollary 2.6. Let he the local time of the standard Brownian motion. Then for all 
A G M, t > 0, and x G M, 

E [exp (A^L^)] = 2 - k|/Vt) + 2$ (|x|/Vt) - 1. 

In particular, 

E[exp(A%^)] <2e^"*/2 + l. 

Proof. Let u{t,x) be a solution to (1.1) with p{u) = \/2 An, n = I, and uo{x) = 1. By the 
Feynman-Kac formula (1.7) with 77 = 2, xi = 0, and X 2 = x. 


E[n(f/2, 0)n(f/2, x)] = E exp(2A^ / 6x{Bl — Bg)ds 


ft/2 


= E 


exp (2A 


rt/2 


^x(.B'Jds 


6 













= E 


exp J (5,(5:)dr jj = E [exp (A^L^)] , 


where B[ is a standard Brownian motion. Then apply (2.11) with A^, t, and u replaced by 
2A^, t/2, and 1, respectively. □ 

For the p-th moments, we have the following bonnd, which simplihes the expression in 
[3]. Denote 

1 if p = 2, 

2 ifp>2. 


Cp 


Theorem 2.7. Let u{t,x) be a solution to (1.1) starting from p G 2 Wh(®)- For all t > 0 
and x G M , the following moment bounds hold: 

(1) If \p{x) I < Lp |a:| for all x G M, then for all p > 2, 

\\uit,x)\\l< // \p\{dzi)\p\{dz 2 )}ClL^{t,x - Zi,X - Z 2 ), 

J J 

where 


zi, Z 2 ) = CpK.*{t,zi,Z2,0]cl^/^Lp)-, 

(2) If |p(x)| > Ip \x\ for a// X G M and if p>0, then 

\\u{t,x)\\l> // p{dzi)p{dz 2 )lC*{t,x - zi,x - Z 2 , 0 -,lp). 


Proof. Part (2) is clear. As for (1), note that the npper bonnds for both 2nd and p-th 
moments share similar forms (see [3, (2.21)]), bnt with different parameters. By the notation 
in [3], Op^o = and Zp = \fpj2. Then replacing the parameter A in /C* by Op^Zp Lp and 
mnltiplying it by a factor Cp, one passes from 2nd moment to p-th moment (see [3, Theorem 
2.4]). □ 

For the alternative proof of Theorem 2.1, we will need the following joint density of the 
standard Brownian motion Bt and its local time at a level a G M, which is by itself 
interesting. When a = 0, it is well known (see, e.g., [5, Exercise 1, on p. 181]) that this law 

is 


P(B;ed'i/,p€dt,)=l^exp 
More generally, we have the following theorem. 


(bl 

2t 


l{^>o}dpdn. 


(2.13) 
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Theorem 2.8. The joint distribution of {Bt, L“) for t > 0 and a G M zs 


P {Bl e dy, L“ e du) = 


\a 

+ \y-a 

+ v 

( (|a 

+ \y-a 


{27lt)^ 

exp 

^ 2t 


f 

e 2 * - 


V^Tlt V 

where sign(a) = 1 if a>0 and —1 if a < 0. 
Corollary 2.9. The law of for t > 0 and a G M zs 


{2a-yY 

e 2 t 


l{sign(a)y<|a|}<5o('y) 


(2,14) 


dz/du, 


P(L‘eAv) = 


TL 

y/wi, 


exp 


{v + |a|)^ 
2t 


+ 


2$ 


|a| 

Vt 


Soiv) l{^>o}dz;. 


Proof. Integrate the right hand of (2.14) for dy over M. For the hrst term, use the integration- 
by-parts formula. For the second term, use the dehnition of <F(-). □ 


2.2 Feynman-Kac formulas for measure-valued initial data 

Let H he a. Hilbert space with inner product (•, •) and let W = {W{h), h G Pf} be a zero 
mean Gaussian process with covariance function K{W{h)W{g)) = {h,g). For any square 
integrable random variable F G let 


n=l 


be its chaos expansion, where /„ G Fr®” (symmetric tensor product) and In denotes the 
multiple stochastic integral. Let L be the generator of the Ornstein- Uhlenbeck semigroup, 
i.e., LF = —J2'^=i''^^nifn)- For any s G M, denote by W’P{H) the completion of if-valued 
polynomial random variables with respect to the norm 


|F| 


I s,p 


= \ \{i-Ly^^F\ 


\LP{n-,H) 


We refer [8] for more details. 

In our framework, H = L^(R+;M") and for any h E H,W{h) = /o°° where 

{Bl, t > 0}, i = 1... ,n, are n independent standard Brownian motions on M. 

Theorem 2.10. For any /i* G Al/f(M), x* G M, hj G if, z = 1,..., n. It holds that 

n 

yiiiW{hi) + Xi) G ©“"’^(M) if a + n/p > n, a > 0 and p > 1. 












Theorem 2.11. Let = {Bl,t > 0}, i ^ I, be one-dimensional Brownian motions belong¬ 
ing to Gaussian space spanned by W, where I is a finite set with m elements. Let L]’^, i E I, 
be the local time of B'^. Suppose that the function f : i—>■ M satisfies that -^f > 0 and 

max sup ||/i + ei,... + em)|| <+ 00 , (2.15) 

for all t > 0, Xi E M., and p >2, where fi = -^f- Then 

f {Ll’^fi ..., G D“’P(M) tfp >1 anda< 1/2. 

In particular, one can choose f{zi, ..., Zm) = exp , A G M. 

Choosing hi = l[o,t], !<*<'«-, Theorem 2.10 implies that 

n 

{fii + i?|) G D““’^(M) if a + n/p > n, a > 0 and p > 1, (2.16) 

i=l 

for uq E Al_f/(1R)- On the other hand, choosing I = {{i,k) : 1 < j < k < n} and B^^T) = 
B^—B^, Theorem 2.11 implies that the random variable Yn dehned in (1.8) belongs to ©"’'^(M) 
if g > 1 and a < 1/2. Therefore, by choosing p close to 1 such that n(l — 1/p) <1/2, and 
choosing q such that l/p+ 1/g = 1, one can see that the Feynman-Kac formula (1.7) is well 
dehned as the duality relationship (1.8), for any uq E AiniM)- This can be proved using the 
fact that Uq * Gi(e, •) converges to uq, as e tends to zero, in the topology of D““’^(R). 

3 Proof of Theorem 2.1 

Proof of Theorem 2.1. The proof consists of the following two steps: 

Step 1. We hrst prove the second moment (2.5). Write Jq in the form of double integrals 

( Jq */C)(f, x) = f ds 

Jo 

By [3, Lemma 5.4] and the notation (2.6), 


dy Kfifi - s,x - y) // p.{dzi)fi{dz 2 )G,,{s,y - zi)G,,{s,y - Z 2 ). 


Gfis, y - zi)Gfis, y- Z 2 ) = Gj./ 2 (s, y - z)Gfi2s, Az). (3.1) 


Denote 


Hit) 


1 A^ pt 

—-|- - 6 $ 

y/Amfi 2v 



(3.2) 
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By Fubini’s theorem and the semigroup property of the heat kernel, 


(Jq */C)(t,x) =A^ // /i(dzi)/i(d22) / ds Gi,{2s,Az) / dy IC{t - s,x - y)Gty/ 2 {s,y - z) 
J Jk2 Jo Jr 

=A^ // yu(dzi)/i(d2:2)Gi//2(^,3; — f) / ds Gi.( 2 s, Az)Fr(t — s). 

J Jr 2 Jo 

Now we will use the Laplace transform to evaluate the ds-integrah By [6, (27) on p. 146], 

exp (—i> 


C[G 2 u{-,x)]{z) = 


— '2 1/2 1x1 y/z) 


2^/Uz 


(3.3) 


By (1) on p. 137 and (5) on p. 176 of [6], 

A^ 1 


C[H]{z) 


+ 


Auz — 2y/17z 


A^ 

2y/i7z {Auz — A^i) 


1 / 1 _ 1 A 1 

2 \2y/i7z — A2 2y/vz + A^y 2y/vz 

A^ ( 1 _ 1 A 

\2-y/Fi — A2 2-01i+A2y 


Hence, 

C[H]{z)C[G2u{-,x)\{z) = fl,-{z,x) - fi,+ {z,x) + f2{z,x) + f3^-{z,x) - f3,+ {z,x), 
where 


fi,±{^,x) 

f2{z,x) 

f3,±{^,x) 


exp (—i/“i/2 |x| ^/z') 
AyJUz {2^/17z ± A2) 
exp (—z/“i/2 |x| y/z) 
Avz 

A^exp (-i/“i/2 |x| ^/z) 
8 uz {2^/uz ± A2) 


Apply [6, (16) on p. 247] with a = z/-i/2|x| and (3 = ±^, 


^ ^[/i,±(-,a;)](f) 


1 J-Aid-l-Alt r 
=—e 21/ ^ iv erfc 

81/ 


\x\ ^ \^Vi\ 

\/Aut a/Iz/ / 


Apply [6, (3) on p. 245] with a = v ^\x\‘^ 


^ ^[f2{-,x)]{t) 



(3.4) 
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Apply [6, (14) on p. 246] with a = u and P = 

Therefore, the ds-integral is equal to 


^ ^ f f f \^\ \ ±l!M+>di r f \x\ . 

± — erfc , — e 2“^ erfc , ± 

Su y \\/4nt / \\/4z/t 


\/4n 


ds G 2 i/(s, Az)H{t — s) = —e 21 . ' + 4^- erfc 


/ |Az| _ A^A 
V\/4z/f \/4i/ / 

1 a 2|A^| A^t f )P\/t \Az\ \ 

—e 2!/ + 4i- $ I-1-L 

2z/ V\/^ V^J' 


Finally, by (3.1), 


Jo{t,x)= // /i(d2;i)/i(d2:i)G'^/2(t,a; - 2;)G2i.(f, A 2 ;). 


This proves the formula (2.5), i.e.. 


\u{t,x)\\l = 


X 


fi{dzi)fi{dz 2 ) G^i 2 {t,x - z) 

. A^ a 2|A2| A^t /A^\/f lA^I 

V 2i/ V 


(3.5) 


Step 2. Now let us consider the two-point correlation function (2.1). Fix t > 0 and 
xi,X 2 E M. Apply (1.6), (3.1) and (3.5), and then integrate over dy using the semigroup 
property: 


X(f,a;i,f,a;2;z/, A) =A^ / dr / dy jju(r,y)jj2G^/2(t - r,x - y)G2u(t - r, Ax) 

Jo Jr 

=A^ dr dy G^/2{t - r,x - y)G2u{t - r,Ax) // y{dzi)iJ,{dzi) 
Jo Jr J Air 2 

/r |A 2 ;| 


where 


„ , -N / A \ a2|Az| A^r /A 

X G^/ 2 {r, y - z) ( G 2 v(r, A^;) + 2 . + 4. $ ( _ 


\/2 


vr 


=A / / y{dzi)y,{dz 2 ) Gy/ 2 {t,x - z) J dr G 2 u{t - r, Ax)H{r, Az), 


A^ _2dM I A* 
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Here, H{t,0) = H{t)] see (3.2). Notice that 


1 _ Afj£| 



t |x| 


-e + "iJ $ I _ 

4z/ \ V 2z/ ^/2Ui) 8z/ 


1 A^lr^l 


= —e 21/ + 41/ erfc — 


t |x| \ 

4z^ \/4i4 y 


Hence, ^ 

H{t,x) = G 2 u{t,x) + 2X^C~^[fi_{-,x)]{t). 

Together with (3.3), we have that. 


C[H{-,x)]{z) = 


exp (—z/ ^/‘^\x\^/z) 
■\f^uz — A^ 


and 




\/AlVz (\/4z/2; — A^) 


By (3.4), 


, ^ / A A \ 1 A^(iAxi+iAzi) A^t / X^\/t |Aa;| + |A^|A 

dr G 2 !y(t — r, Ax)H{r, Az) = —e 21/ + 41/ $ I - - - -- . 

V v2za Aj^vt J 


2v 


Therefore, 

¥.[u{t,Xi)u{t,X2)] = Jo(t,Xi)Jo(t,X2) 


/i(dzi)/i(dz2) Gi// 2(^, a; - 2r) 


A^ A^(|Ax| + |Az|) 

X —e 21/ 

2v 



Ax + A^; 


y/Aui 


This completes the proof of Theorem 2.1. □ 

4 Proof of Theorem 2.8 

Proof of Theorem 2.8. The case when a = 0 is covered in (2.13). Let g : M. x M+ M be a 
smooth and bounded function. Denote the joint density of {Bt,L^) by fa,t{,x,y). Then 

poo p 

E[^(Bi,^“)]= / d|/ / dx g{x,y)fa,t{x,y). 

Jo Jr 

First assume that a > 0. By the reflection principle (see [10, p. 107]), the density of 
Ta := inf (s > 0 , = a) is 

fxa (s) = exp ( j ’ for s > 0 and a > 0. 


\/2 


TTS"' 
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Let be the standard Brownian filtration. By the strong Markov property, 

E |j(B„ L?)l = E [E (9(B., Lr)\FT.} 1|t.<<)] + E [E (g(B„ Lt)\TT.) 1|t.><)] 

= / ds/T.(s)-B 9 (Bt.-, + a,L'l_,) + E (g(B,, 0 )l(T.>t)) 

Jo 

=■■ h + h, 

where I3t is a standard Brownian motion and is its local time at the level 0. We hrst 
compute Ji- By (2.13), 


h = 


dydv g{y + a,v) / ds 


V2 


TTS'^ 


: exp — — 


\y\+v 


j j ^ , ^a+\y\+v 

dydv g{y +a, v) - exp 


2s / ^27r(t — s)3 
a+ \y\+vf 


exp 


+ vY 


2{t - s) 


2t 


where in the last step, we have used the fact that the densities form a convolution 
semigroup, namely, = fT^+b (see [10, p.l07] or Lemma 4.1 below for a short proof). 

As for I 2 , let Mt = sup^<ii?s. Then, by the joint law of (see, e.g., [10, Ex. 3.14 

p. 110]), 


/■“ r™ \/2(2m - v) 

l 2 = E[g{Bt,0)l{Mt<a}] = dm dy g{y,0) - ’ exp 

JO J—00 V 


dyg{y,0)- 


'y+ 

1 




(2m — yY 
2t 

(2m — yY 


2t 


c 


(2y^-y)^ (2a-y)^ 

e 2t — e 2t 


\/27ri V 

where y+ = inax(y, 0). Note that 2y+ — y = \y\- Combining these two parts, we see that 


fa{y,v) = 


a +\y — a\+ V 


{27itY 

1 ( 


■ exp 


{a+\y - a\+ vf 
2f 


\/2nt V 


_y_ _(^a-yj_ 1 r / \ 

e - e 2 t l|j^<„|(5o(n) 


4i;>0}- 


If a < 0, by symmetry, let Bt = —Bt and be its local time. Clearly, the density of 
{Bt)L^°‘) is Therefore, the law of {Bt,L‘l) is f\a\{—y,v). Then use the fact that 

I — 1/ — |a| I = \y — a\ and |2|a| + y\ = \2a — y\. This completes the proof of Theorem 2.8. □ 


Lemma 4.1. For t > 0 and a,b Y d, the following integral is true 


\ab\{s{t-s)) ^/^exp^^-^ 


62 


2{t — s) 


ds = 


|a| + |6| 

Y27iF 


■ exp 


i\a\ + \b\f 

2t 
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Proof. Denote the integral by I{t) and let gait) = |a| t 3/2ga^/(2d_ By [5^ (^28) on p. 146], 

C[ga]{z) = V^e-^\^\^P (4.1) 

So, C[I]iz) = C [ga] iz)C [g}f\ {z) = Then use (4.1) for the inversion. □ 


5 An alternative proof of Theorem 2.1 


If the initial data fi is such that fiidx) = uoix)dx, where uq is a bounded measurable function, 
then we can use the Feynman-Kac representation (1.7) with n = 2 and Theorem 2.8 to give 
an alternative proof of Theorem 2.1. 

An alternative proof of Theorem 2.1. For initial data specihed above, E (-u(t, xi)uit, X 2 )) ad¬ 
mits the Feynman-Kac representation (1.7). Let 

Wl,-.= 31,-31, and Wl,-.= 3l, + 3l,. 

Note that Wf and Wf are two independent standard Brownian motions. Then, 


E (M(f, Xi)uit, X2)) = E {xi + 2"^[PF2^t + w^^,]) Uq (xa + 2.-\Wi^, - W^^,]) 

X exp 

= [ dz G 2 ^it, z)eI Uq {xi + 2~^[Wf^, + z]) Uq {X 2 + 2~^[z - Wf^,]) 

Jr y 

X exp ds 

Notice that 
exp 


= exp (hF,^) ds^ = exp 


2l't I 5 


where L“ be the local time of Wf at the level a. Then the expectation in the above integrand 
becomes 


/A' 


E(tio(2:i+2 ‘|*r2„i + 2|)!io(i:2 + 2 ‘|2- W^2j)exp ( —L£, " 


(5.1) 


By Theorem 2.8, this expectation is equal to 


, , z + y\ ( z — y 

dy Mo a;i H -) mo I + 
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^ \y — Axl + \ Ax\ + V 
X / df - , -exp 


(ll/— Aa;| + f + |Aa;|) ^ 


(5.2) 


'0 


4 \/^ \ 4z/t 2u 

dy uo{xi + {z + y)/2)uo{x2 + {z - y)/2) [G 2 y{t, y) - G 2 u{t, 2Ax - y)]. 


J sign(Ax)y<| Aa:| 

By integration by parts, the dt;-integration in (5.2) is equal to 


1 \^v (ly-Axl + IAxId-u)^ 

g 2i^ g 4i/t 


\/ A-KUt 


V —0 \ 2 POO 1 o 9 

A / 1 X^v (|y-Ax| + |Ax|+-u)^ 

+ — / dw , e 2^ e 


V = 00 2:1^ Jq 


\/ Ain't 


, , , , , ,, X^t A2(|y-A:.| + |Ax|) / „ HT lu — Axl + lAxI 

=G,„((, I, - Ax| + |Ax|) + - 4 (^A^- - 

Denote the d|/-integral in (5.2) by I. By symmetry, we will only consider the case where 
Ax > 0. In this case, the d|/-integral is from —oo to Ax. By the change of variables 
y' = 2Ax - y, 


'*Ax 


dy Uo{xi + {z + y)/2)uo{x2 + {z - y)/ 2 )G 2 y{t, 2Ax - y) 

POO 

= / dj/' + (z + v')/2)uo{x2 + (z- !/')/2)G2„(i, y'). 


(5.3) 


' Ax 


Hence, 


1= dyuo{xi +{z + y)/2)uo{x2 + {z-y)/2) [G 2i/(^) y) l{2y<Ax} 2/)1 {j/>Ax}] 


Notice that 


G 2 v{t, y)l{y<Ax} - G 2 u{t, y)^{y>Ax} + G 2 u{t, \y - Ax\ + Ax) 

= [G 2 ,y{t, y) + G 2 u{t, 2Ax — y)] l{y<Ax}- 

Therefore, by (5.3), we know that 

I dy Uo{xi + {z + y)/ 2 )uo{x 2 + {z - y)/2) 

J'R 

X (G 2 y{t, y)l{y<ixx} - G 2 u{t, y)l{y>Ax} + G 2 u{t, \y - Ax\ + Ax)) 

= / dyuo{xi + {z + y)/2)uo{x2 + {z-y)/2)G2u{t,y). 

Jr 

Combining these calculations, we have that 

E {u{t, xi)u{t, X 2 )) = dzdy Uo{xi + {z + y)/2)uo{x2 + {z - y)/2)G2u{t, z) 
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X 


, , A^(|y-Ax| + |Ax|) 

G2.(t,2/) + —e- 2. $ AW-^- 



2 / ^ b-Ax| + |Ax| 


\/2i/t 


Finally, by change of variables 2:1 = Xi + ( 2 ; + 1/)/2 and ^2 = X 2 + ( 2 ; — 2 /)/ 2 , we have that 
E {u{t, xi)u{t, X 2 )) = dzidz 2 uo{zi)uo{z2)2G2u{t, {zi + Z 2 ) - (xi + 0 : 2 )) 


X G 2 u (t, Az - Ax) 


f XH A^l^i - 2 : 2 ! + |Aa;|A / 2 |^i - 2^21 + |Aa;| 


2.""P Id. 


2 . 


$ A^ 



2 . 


\/2.t 


Notice that 2G2u{t, (zi + Z 2 ) — (xi + X 2 )) = Gu/ 2 {t,x — z). This proves (2.1) (see (2.8) and 

( 2 . 10 )). □ 


6 Feynman-Kac formula for measure-valued initial data 

6.1 Initial data part (Proof of Theorem 2.10) 

We need some lemmas. 

Lemma 6.1. For any t > 0, s > 0, x G M, ?/* G M, i = 1,... it holds that 


Gi{s,x + z - yj)dz < (p+lfG. —^—eii+q nGi((p+ 

fi Vps + t 


Proof. Denote y = {yi + ■ ■ ■ + yp)/p. Notice that 

JJ Gi{t, z - yj) = (27rt)"^/^ exp (- j exp 


i=i 


2t/p 


= {2TTt)-^^-^^/^p-^/^G,{t/p,z-y) 


yl + '-' + yl-py"^ 

2t 


exp 


yl + '-' + yl-py"^ 

2t 


Hence, by the semigroup property of the heat kernel, the dz integral is equal to 

p 


/ dzGi{s,x + z) n Gi{t,z-yj) 

Ar 

=(27rf)“^^“^^'^^p“^'^^Gi(s + t/p, X + y) exp 


yl + '-' + yl-py"^ 

2t 
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= {27r)-P/H-^P-^^/^{ps + exp (^ 

<{27r)-P/H-^P-^^/^{ps +exp 
=(27r)“P/^t“(P“^^/^(ps + t)-V2g2(f;^ exp 

exp 


p{x + yf yl + --- + yl-p y" 


2{ps + 1) 2t 

pf - 2 px^ yl^ - \-yl-py 


A{ps +1) 2t 

pf{t + 2ps) yl + --- + y.^ 


At{ps + 1) 


2t 


4)fe? + - + !/S) 


<(27r) + exp 


4t(ps + 1) 2t 

yl + --- + yl 


2{p+l)t 


-{p + l)^/^e2(p»+‘) n Gi((p + l)t,2/i). 


2=1 


ps +1 

This proves Lemma 6.1. 

Let Tt be the Ornstein Uhlenbeck semigroup of operators on L‘^{Q) with generator L. 
Lemma 6.2. For any pi G Aini^), hi G H, and Xj G M, i = 1,..., n, it holds that 

/ n \ n 

Ttil[p^ {xi + W,{hi)) = J] [p, * Gi(|h,|2(l - e-^O, •)] + e-^Wi{hi)), 


□ 


<2=1 




2=1 


for t > 0, where Wi are i.i.d. zero mean Gaussian processes Wi = {Wi{h),h G H} with 
covariance function E,(Wi{h)Wi{g)) = {h,g). 

Proof. Fix e > 0. Let /ij,e(x) = (/ij * Gi(e, •))(3^)- By Holder’s inequality, 


E 


Xi 


< 


2 = 1 


_2 = 1 

Notice that for all p > 1, 

¥.[\pi^^{Wi{hi) + Xi)f] = [ dz Gi{\hi\'^,x + z 


l[E[pi^,{W,{hi)+x,f^] 


2nll/« 


Gi{e, z - y)pi{dy) 


< / \Pi\{dyi).. .\pi\{dyp) / dz Gi{\hi\^,x + z)YlGi{e,z - yj). 
Jm.p Jr 

By Lemma 6.1 with t = e and s = 

px^ 

E [\p^,,{WM + Xi)\P] < (p + l)PAeAFWTo [(|^.| ,, G'i(e(p + 1), •))( 0 )f , 
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which is hnite because //* G Hence, nr=i+ Xj) G -h^(H), and one can 

apply Mehler’s formula (see, e.g., [8, Section 1.4.1] to obtain that 


Tt[\{^^e{x, + W,{hi)) 


v2=l 


/ 


/ii,e + e ^Wi{hi) + \/l - W-{hi)^ 

n 

Wliie {xi + e~^Wi{hi) + y^) Gi (|hi|^(l - e~^^),yi) dyi 
i=l 


n * Gi{\hi\\l - e-2*), •)] {X + e-^W,{h,)) 

i=l 


= ll[f^i* G,{e + \h,\\l - e-^0, •)] (x. + e-^W,{h,)). (6.1) 

i=l 

Finally, Tt (nr=i + Wi{hi))) is the L2(H)-limit of Tt (nr=i and this 

limit can be obtained by sending e to zero in (6.1). This proves Lemma 6.2. □ 


Proof of Theorem 2.10. Without loss of generality, assume that /Xj > 0. Following [9], we 
write (/ — L)““/^ in the following form 

poo 

(J - L)-“/2 = F(a/2)-i / e-H^/^-^Ttdt. 

Jo 


By Lemma 6.2, 


i=l 

poo ^ 

F(a/2)-i / e-‘r/2-i JJ * G'i(|h,|^(l - e"^*), •)] (e-‘W,(h,) + xO)df 


2=1 


< F(a/2)-^ / e-HG'^-^ 

Jo 


n h *G,{\hJ\l- e- 2 *),.)] (e-*W,(h,) + X,)) 


i=l 


P 

df. 


( 6 . 2 ) 


Now 


n [p, * G'i(|h,| 2 (l - e- 2 *),.)] (e-*W(h,) + X,)) 

2 = 1 


=E 



(xj + e ^Wi(hi)-y)'^ 

2|hiP(l-e^t) 


>/27r|hi|2(l - e-2*) 


/^i(d2/) 



p 


p 
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where 


„ / n \ p p n 

/ dzi...dzn[Y\Gi{Si,Xi +Zi) L nn Gi{Ti,Zi - yij)ni{dyij) 

Jr- Vi=i / j=i i=i 

'' p n \ ^ r ^ 

nn n / + ^i) n “ Vij) 

u=l i=l / i=l 7=1 


^i = |hjp(l —e and Si = e 
By Lemma 6.1, the dz* integral is bounded by 


dz, G^{Si, Xi + n Gi(T,, - y,p <{p + J] 


..2 P 


J = 1 


<T. 


(l-p)/2 


+ T, 

|hdP 


i=i 


,2 p 


-7=^=^e JJ Gi (I hi 12 (p + 1), i/i ■), 

Vp^i + Ti 


where we have applied the inequality 


G^l((P + 1)^7) Vij) < 


|^Gi((p + l)|hi|^, yij). 


Hence, 


nh*Gi(ih,|2(i 

i=l 


e-2*),-)] {e-^W,{h{) + X,)) 


< 


T„. 


|^^|g 2 (pSi+Tp 


n (p^^+T,)V(2p) 


'^o((p+l)|hi 1^, 0). 


By substituting the above upper bound into (6.2), we see that the integral in (6.2) converges 
provided that 


a -I _|_ _7T^ n 

t 2 ■^2p 2 < oo^ 


0+ 


where we have used the fact that 


pSi + Ti>\hi\^pe and Ti = 2\hi\^{t + 0{t^)). 

Therefore, a + n/p > n. This completes the proof of Lemma 2.10. □ 

6.2 Local time part (Proof of Theorem 2.11) 

Proof of Theorem 2.11. This is a slight extension of [1, Theorem 1]. The proof consists three 
steps. Let be the local time of the standard one-dimensional Brownian motion. 
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Step 1. Fix e > 0. Denote 


Define 


1 if y > X + e, 

F,^:,{y)= {{y-x + e)/{2e) ii\y-x\<e, 

0 ii y < X — e. 

N^{x,t) = [ F^^^{Bs)dBs and N{x,t) = [ l{Bs>x}dBs 


By Tanaka’s formula, 

L^ = (5t-x)+-(-x)+-iV(x,f). 

Notice that L{x,t) G [0,t] is a bounded random variable. On the other hand, let 

Fit ■= {Bt - x)+ - (-x)+ - iV,(x, t). 

Notice that 

N[x + e, t) < iVe(x, t) < N{x — e, t). 

Using the fact that for all x G M and y > 0 

x"*" < (x + y)~^ < x~^ + y and x"*" > (x — y)~^ > — y, 

we see that 

< (^Bf — (x + e) + e)~*~ — (—(x + e) + e)~*~ — iV(x + e, t) < + e, 

and 

Lit > {Bt — (x — e) — e)"'' — (—(x — e) — e)'^ — N{x — e,t) > — e. 

Hence, it holds that 


Lf-" - e < Lf i + 


e. 


(6,3) 


Step 2. Now assume that e G (—1,1). By telescoping sum, 

- f{z[,...,z’^)\ < \f{zi,Z2,. . . ,Zm) - f {z[, Z2, ■ ■ ■ , Zm)\ 

+ \f{z[,Z2,Z3, ...,Zm)- f{z[,Z2,Z3, . . . , Zm)\ 

+ ... 

+ ..., <-i, Zm) - f{z[,..., zl_^, z'J I . (6.4) 

Because > 0, the hrst term to the right of (6.4) satishes that 


\f{Zi,Z2,...,Zm) - f{z[,Z2,...,Zm)\ < (|/l (^1, ^2, • • • , ^m) | + \fl{z[, Z2, . . . , Zm)\) 1^1 - z[ 


II- 
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Now replace Zi and z[ by U := and L\ := i = 1,... ,m, respectively. Denote the 
quantity in (2.15) by 

Cp . Cp(t^ Xi,..., Xm). 

Let 1/r + l/g = 1, q > 2. By Holder’s inequality, the hrst term on the right hand side of 
(6.4) satishes that 

||/(L^L^...,L-)-/(4^L^...,L-)|■ 


p 

<2 sup + ei, ,..., L'^)\ 

o n 11 r^i r*i 11 

— 11-^4 ^e,t\\pr ) 


I I 


\pq 


\pq 

1^1,3:1 _ ^1>®1 I 


\pr 


\pr 


where we have used the fact (6.3). This inequality is true for all the n terms on the right 
hand side of (6.4) under the same replacements. Therefore, 

m 

11 /(L‘,,,,, L”) - /(l; ,,.., ir) 11 < 2 C„ ^ 11 if* - L311 < 2mC„£■/^ (6,5) 


i=l 


where the last inequality is due to (2.8) of [1]. 

By Holder’s inequality with l/r + l/g = l, we see that 


m 

E [\\Df{Ll ,.... L”)11^] < 2t-‘ ^E (|/,(P..... Lr)r | |bl:ji i;) 

i=l 

m , 

<2 t-‘J]E(|/,(L;,.,,,Lr)r)‘''''E(||BL:j||^) 

i=\ 

m 

£2’’-‘CJ,^||DL3 


i,a:| IP 

5^ \ \lp^{Q-,H) ’ 


2=1 


where we have applied (6.3) and (2.15) as before. Hence, by (2.9) of [1], 

m 

2=1 

< 2 Cpq V for all g' > 1. (6.6) 

Therefore, by the same arguments as the proof of Theorem 1 in [1], we see that (6.5) and 
(6.6) imply that /(T^,..., L™) G D"’P(M) for all p > 1 and a < 1/2. 


Step 3. In this hnal step, we need to verify that f{zi ,..., Zm) = exp ) satishes 

the condition (2.15). By Holder’s inequality, we have that, for all ej G (—1,1), 


+ ^1? • • • 5 L — A 




exp ( + ej 

1=1 
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j=i 


mp 


< JJ||exp (a^L: 

i=i 


2 rj,^j+<^j 
t 


mp 


By Corollary 2.6, 


exp ( 


= E 


mp 


exp ( 


< 


2gA'‘m2p2t/2 ^ 


1 

mp 


Therefore, 


W £ / T l.xi+ei I T n 

max sup \\ji[W +ei,...,Li 

,e(-i.i 


m,a^rn+en 


l<i<m^^,g(_l 1) 


+ e^)||^ < AV^™ 


2gA4m2p2t/2 ^ ^ 


This completes the proof of Theorem 2.11. 


□ 
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